Anomalies in electrostatic calibrations for the 
measurement of the Casimir force in a sphere-plane geometry 
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We have performed precision electrostatic calibrations in the sphere-plane geometry, and observed 
anomalous behavior. Namely, the scaling exponent of the electrostatic signal with distance was found 
to be smaller than expected on the basis of the pure Coulombian contribution, and the residual 
potential found to be distance dependent. We argue that these findings affect the accuracy of the 
electrostatic calibrations and invite reanalysis of previous determinations of the Casimir force. 
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Over the last decades the Casimir force [l[ has met 
increasing popularity as a macroscopic manifestation of 
quantum vacuum with its relevance spanning from 
nanotechnology to cosmology 0]. With the claimed ac- 
curacy of recent experiments ranging from 15% in the 
arallel plane case [J| to 0.1-5% in the sphere-plane case 
S S 13 1 it also provides constraints on the existence of 
forces superimposed to the Newtonian gravitational force 
and expected in various unification attempts [lo| . 

Concern has been raised that previous analyses in the 
measurement of the Casimir force have overlooked pos- 
sible influence of residual electric effects (the so-called 
patch effects) , which could mimic the Casimir force [llj . 
Residual electric effects are known to play an impor- 
tant role in the measurement of van der Waals force be- 
tween macroscopic bodies, where corrections based on a 
model for work function anisotropics and their associated 
patch charges have been discussed [121 ■ Some of the early 
studies focusing on various adhesion and friction surface 
forces can be found in jl3, 14, 15j. More recently, Stipe et 
al. argue that the presence of an inhomogeneous tip- 
sample electric field in an AFM type experiment is dif- 
ficult to avoid, thereby imposing a significant limitation 
to the accuracy of force measurements. This issue indi- 
cates a common problem in electrostatic measurements 
over a range of different experimental circumstances and 
distance scalings, and calls for more attention on calibra- 
tion procedures, as the accuracy claimed in Casimir force 
measurements inherently depends on the quality of the 
corresponding electrostatic calibrations. 

We report here an electrostatic calibration procedure 
for a sphere-plane geometry in which the electric signal 
is studied at all explored distances. Our results are ob- 
tained in a range of parameters that interpolates between 
the two previous sets of sphere-plane measurements. We 
use a spherical lens with a large radius of curvature, simi- 
lar to the experiment performed in , while at the same 
time exploring distances down to few tens of nanome- 
ters from the point of contact between the sphere and 
the plane, similar to more microscopic setups using mi- 



croresonators 0, 0, Q • The measurements reveal anoma- 
lous behavior not reported to date. 

Our experimental setup is an upgrade of the previ- 
ous arrangement for measuring the Casimir force in the 
cylinder-plane configuration [17| . Force gradients be- 
tween a spherical lens and a silicon cantilever are de- 
tected by measuring the shift in the mechanical oscil- 
lation frequency of the cantilever. A schematic of the 
apparatus is shown in Fig. [TJ The cantilever is electri- 
cally isolated and thermally stabilized by a Peltier cooler 
to within 50 mK. The motion of the cantilever is moni- 
tored by a fiber optic interferometer [3] positioned a few 
tens of microns above the cantilever. The output signal 
from the interferometer is put through a single reference 
mode lock-in amplifier and is fed back into the piezo- 
electric actuator driving the cantilever motion, forming 
a phase-locked loop (PLL) with an optimized phase angle 
around 30-40 degrees The measured frequency at a 
typical vacuum pressure of 1.6 x 10~^ Torr is consistent 
with the predicted frequency of the fundamental flexural 
mode of the cantilever, around 894 Hz. The stiffness k 
of the resonator is estimated to be 5.4 x 10"^ N/m, 10^ to 
10^ times higher than that of typical cantilevers used in 
atomic force microscopy. Our cantilever flexes at most 1 
nm, allowing very small gaps to be probed, at the cost of 
an overall lower force sensitivity that limits the maximum 
explorable distance. No compensating external voltages 
are needed to prevent snapping. 

In general, the square of the measured frequency Vm 
of a cantilever under the influence of generic forces is 

ylid, V) = ylid ^ oo) - l^ylid, V) - Az.,2(d), (1) 

where Vp is the cantilever's natural flexural frequency, 
Az/g is the frequency shift due to externally applied volt- 
ages V at different gap separations d, and Ai^^ is the fre- 
quency shift subject to distance-dependent forces of non- 
electrostatic nature, for instance the Casimir force. For 
the sphere-plane configuration and for our choice of pa- 
rameters, the proximity force approximation (PFA) [20| 
holds, the electric force gradient is F^'j = Tiei^RV^ / d"^ , 



FIG. 1: (Color online) Schematic of the experimental setup. 
The silicon resonator, of length L = (22.56 ±0.01) mm, width 
w = (9.93 ± 0.01) mm, and thickness t — (330 ± 10) fim with 
the physical mass rup = (1.72±0.05) X 10"" kg, is opposed by 
a spherical mirror with radius of curvature R = (30.9 ± 0.15) 
mm and diameter a — (8.00 ± 0.25) mm. The mirror is 
mounted on an aluminum frame connected to two motorized 
actuators allowing for coarse translational motion, plus an ad- 
ditional piezoelectric transducer for fine translational motion 
driven with a bias VpzT- Both the spherical mirror and the 
surface of the cantilever are coated by thermal evaporation 
with a 200 nm layer of gold with 1-2 nm rms roughness. 



and Eq. ^ can be parameterized in the presence of a 
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contact potential Vc as 



kei{V - V,y, where 



1^0 = t^p — , a parabola whose maximum is reached 
when the applied voltage equals to Vc- The parabola 
curvature kd = eoR/ATrrrioscP reflects the cantilever re- 
sponse to externally applied electric forces at a given dis- 
tance. All three parameters, iyQ,kci, and Vc are simulta- 
neously evaluated for every step in a data sequence and 
so they can be plotted as a function of gap separation. In 
what follows, we scrutinize these parameters and discuss 
the systematic issues that must be addressed to validate 
a consistent analysis of the Casimir force. The value of 
fcei as a function of distance directly characterizes the 
system's response to an applied bias and sets the basis 
for residual distance-dependent force analysis. The gap 
distance cannot be known with sufficient accurately prior 
to a force measurement [lol |. In a typical experiment, the 
gap is varied by the voltage applied to the PZT (Vpzt) 
and consequently feci is a function of relative distance (i.e. 
of the applied Vpzx)- This requires an additional fitting 
parameter Vp-^-j , which would cause contact and must be 
inferred from fitting the function 



fcci(VpzT) = a{Vp2T - Vpzt)" 



(2) 



where a is a calibration factor containing the effective 
mass of the cantilever through a = eoR/ATrniesP^ , and 
(3 is the conversion factor translating the PZT voltage 
to the actual distance unit. This is common practice 
in most of the Casimir force measurements in which the 
sensitivity of the apparatus and the zero distance are ex- 
tracted from an electrostatic calibration. Depending on 
apparatus type, extracted quantities can be cither spring 



FIG. 2: (Color online) Curvature coefficients data k^i versus 
applied voltage Vpzx and best fit with Eq. ([5}, but leaving the 
exponent as a free parameter for four different experimental 
runs. Given 4% uncertainty in kei , the reduced bltb 1.0, 
0.8, 1.2, and 7.0 with the exponent left as a free parameter 
for Run 1, Run 2, Run 3, Run 4, respectively. Fixing the 
exponent to be -2 increases these values to be 16, 7.7, 6.9, 
and 37. The data follows a pure power law in all explored 
distances. The absolute distances can be assessed from the 
asymptotes of the curves V^pztj which are, in the progressive 
order of run, 43.12 ± 0.01 V, 56.77 ± 0.02 V, 56.99 ± 0.02 V, 
and 56.66 ± 0.02 V. The different values of Vp^T in the four 
runs, taken in different weeks, are due to the rearrangement 
of the sphere location obtained using the actuator drivers. 



constant @ , torsion constant 0, 0, Si , or effective mass 
like in this case [i, 17|. Note that the absolute distance 



can be expressed in two ways: d(VpzT) = /^(VpzT" Vpzx) 

1/2 

or d{VpzT) = f3[oi/kci{VpzT)] ■ Therefore, the absolute 
distance can be inferred either from the asymptotic limit 
^PZT ^^'^ function or from the calibration factor a 
of the same function, indicating an interdepcndency of 
the two physical parameters appearing in Eq. 

The above procedure and all subsequent analysis is in- 
applicable if the data fail to follow the inverse square law 
ofEq. Surprisingly, our experimental data from four 
separate sequences follow a power law similar to Eq. 
but with exponents -1.70 ± 0.01, -1.77 ± 0.02, -1.80 ± 
0.01, -1.54±0.02, far from the expected value of -2. Fig- 
urc[2]shows plots of k function of VpzT , with the ex- 
ponent left as a free parameter. With the apparatus sta- 
tionary, we observe a random measurement uncertainty 
of 4% in kei superimposed on a long time scale (compa- 
rable to the duration of one data run) drift, which may 
be due to a systematic, thermally induced change in gap 
size up to ±200 nm. Given the 4% error, the reduced 
is near one when the exponent is fit in the -1.7 to - 
1.8 range compared to around 10 for fixed -2.0 exponent. 
Adherence of our data to a strict power law from the far- 
thest to the closest approach in all four cases implies the 
drift was not severe during our data runs, except possi- 
bly in the case of Run 4. To check the stability of the 
electrostatic result ruled by the unexpected power law, 
the following test has been conducted. We repeat fitting 
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FIG. 3: (Color online) Stability test of fit parameters in 
Eq. (I2J for Run 1 with the exponent taken to be -1.7 (black, 
circles) and -2 (red, diamonds). Both calibration factor a 
(top plot) with units of Hz^ m^'^ for the -1.7 exponent, 
and Hz^ m^ for the -2.0 exponent, and the zero distance 
bias Vpzx (bottom plot) remain stable within 3% in the flat 
region at the closest separation for the exponent -1.7, while a 
large spreading in the two parameters is evident for the expo- 
nent -2. Note also the correlation between the fit parameters 
reflecting their interdependency as discussed in the text. 



the data, fed versus VpzT, starting with few points at 
the largest distances and by progressively including the 
data point corresponding to the closer distances. In this 
way the intrinsic instabiUty of the fit will manifest itself 
through a systematic variation in the fit parameters sub- 
ject to a choice of data points included. Fig. [3] shows the 
result of the stability test for Run 1 . The electrostatic fit 
with the exponent of -1.7 displays a stable region at the 
shortest separation gap corresponding to the inclusion of 
most data points. Excluding data points simply makes 
the error bar larger, while the fit parameters remain sta- 
ble. Forcing the exponent to be -2, on the other hand, 
leads to unsatisfactory calibration results, as the param- 
eters start to drift immediately from the outset and do 
not display any stable region. More than 40% variations 
in both the calibration factor and the zero distance are 
found in this case. The instability is already present at 
the largest distances, indicating that the source of devia- 
tion from pure Coulombian behavior, for instance due to 
patch potentials, is still effective. Runs 2, 3 and 4 exhibit 
similar behavior. In order to further test the consistency 
of the electrostatic calibration, we have compared the ab- 
solute distances obtained from the two methods outlined 
above, getting a mutual agreement within 2 % for the fit 
with the -1.7 exponent at all distances, and deviations of 
20 % if the fitting exponent is instead fixed to be -2. 

The unexpected power law poses a significant limit 
on the validity of our electrostatic calibration. Here, 
we briefly examine some hypotheses which could potcn- 
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FIG. 4: Dependence of the contact potential versus sphere- 
plane separation. Two noticeable behaviors are a linear de- 
pendence on the separation (Run 1 and Run 2) and a settle- 
ment to a constant contact potential at the smallest explored 
distances (Run 3 and Run 4). In Run 1 the uppermost point 
of the sphere is different from the other, successive three runs 
as the latter were taken after having tilted the sphere by 0.05 
radians, to check the sensitivity to local surface details. Note 
that in Run 1, Vc is positive while it is negative and converges 
to Vc — -150 mV at the smallest separations for other runs. 



tially explain a deviation from the expected power law. 
Static deflection of cantilever: The spring constant of our 
cantilever is extremely stiff (about 5400 N/m). Using 
Hooke's law, a deflection experienced by the cantilever 
due to an electrostatic force at 100 nm with an applied 
voltage of 100 mV is less than 0.2 A. Hence, the static 
deflection should play little role. Thermal drift: Even 
though the temperature of the cantilever is actively sta- 
bilized by a Peltier cooler, the rest of the system is still 
subject to global thermal variation. In order to see this, 
we have measured fcei with respect to time at a nomi- 
nally fixed distance. In the worst circumstance, the gap 
separation during the course of measurements can drift 
as much as 200 nm in either direction. Although such 
a drift could in principle affect the inferred exponent, a 
highly unlikely non-linear monotonic drift would be nec- 
essary to account for the consistently observed anomaly 
in each independent run. Non-linearity of the PZT trans- 
lation: The linearity of the PZT translation has been 
tested under a number of different circumstances. Notice 
that the translation intervals between the data points in 
each of the runs shown in Fig. [2] are completely random. 
Yet, all of the runs obey a specific power law in all dis- 
tances. The PZT was also independently calibrated by 
fiber optic interferometer with a consistent, linear actua- 
tion coefficient factor (3 = 87±2 nm/V. Non-linear oscil- 
lation of cantilever: The cantilever is driven at resonance 
in a phase-locked loop, a routine technique adopted by 
many groups 
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[201 ■ Higher order terms in the 
force expansion should produce higher harmonics of the 
drive frequency. Then the assumption that the frequency 
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shift is simply proportional to the gradient of the exter- 
nal force F'{d) could lead to erroneous assessments of kei, 
eventually affecting the exponent. We have not observed 
higher harmonics in the frequency spectrum of the res- 
onator. Surface roughness: The deviation from geomet- 
rical ideality and its influence on the local capacitances 
[23! could in principle play a role especially at the small- 
est distances. However with the measured rms values 
of roughness for the two surfaces we find the corrections 
negligible, as we discuss in detail elsewhere [24]. 

After discarding the hypotheses above, we may con- 
sider the effect of the patch surface potentials. Patch ef- 
fects are expected to induce deviations from the Coulom- 
bian scaling with distance and we have found an- 
other anomaly corroborating this hypothesis. Indeed, 
our technique for obtaining the parabola curvatures at all 
distances reveals that the contact potential depends on 
distance, as shown in Fig. |4l In runs 1 and 2, the contact 
potentials appear to be distance-dependent. This finding 
calls for a more careful analysis of previous experiments 
in which the determinations of the contact potential have 
been performed at relatively large distances, typically 
above 1 fj,m. To see this more clearly, we evaluate the 
equivalent voltage Vcq = [ir / d) ^/hc/^QtQ necessary to 
mimic the Casimir force at a given distance, as first dis- 
cussed in |25| . At 1 /xm, the magnitude of the Casimir 
force is equivalent to that from an uncompensated volt- 
age of 10 mV between the two surfaces, which should be 
compared to variations of 90 mV and 50 mV of the mea- 
sured contact potentials in Run 1 and Run 2. respectively. 
The determination of the contact potential at large dis- 
tances, as usually performed in various experiments, may 
therefore lead to a spurious signal of electrostatic origin 
if the contact potential at smaller distances is only partly 
compensated by an external countcrbias. 

Once the anomalous scaling exponent and the distance 
dependent residual voltage are taken into account, one 
can look for the distance dependence of i/q. A distant- 
dependent frequency should signal forces such as the 
Casimir force and/or forces due to surface potentials re- 
lated with patch effects [Ol ■ Disentangling forces of dif- 
ferent origins becomes exceedingly complex with the un- 
expected power law found in the electrostatic analysis. 
Preliminary residual fitting indicate that the exponents 
obtained for Run 1 and Run 3 arc —2.12 ± 0.16 and 
—3.64 ± 0.25, respectively, systematically smaller than 
-4 expected for the Casimir force (23 |. 

Finally, we emphasize that the study of the Casimir 
force should be regarded as an extension of previous van 
der Waals force measurements with AFM techniques, 
since the underlying physics governing the short range 
forces between closely spaced bodies should be the same. 
Some of the systematic effects discussed here have been 
extensively discussed in the AFM literature [22, 2^, 27 1. 
Although we cannot draw substantial conclusions about 
the observation of the Casimir force itself from our data, 



the experimental procedure outlined in this report should 
provide an optimal strategy to handle contact potentials 
at all distances. Apart from looking for confirmation of 
the observed anomalies in other experimental setups, we 
believe that our findings call for a reanalysis of previous 
Casimir force experiments in the sphere-plane geometry, 
a check of the claimed accuracy, and for the systematic 
control of the effect of uncompensated patch potentials. 
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sity, 217 Prospect Street, New Haven, CT 06520-8120. 
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